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We present reduced-dimensional stochastic projected Gross-Pitaevskii equations describing regimes of con¬ 
finement and temperature where a ID or 2D superfiuid is immersed in a 3D thermal cloud. The projection 
formalism provides both a formally rigorous and physically natural way to effect the dimensional reduction. 
The 3D form of the number-damping (growth) terms is unchanged by the dimensional reduction. Projection 
of the energy-damping (scattering) terms leads to modified stochastic equations of motion describing energy 
exchange with the thermal reservoir. The regime of validity of the dimensional reduction is investigated via 
variational analysis. Paying particular attention to ID, we validate our variational treatment by comparing nu¬ 
merical simulations of a trapped oblate system in 3D with the ID theory, and establish a consistent choice of 
cutoff for the ID theory. We briefly discuss the scenario involving two-components with different degeneracy, 
suggesting that a wider regime of validity exists for systems in contact with a buffer-gas reservoir. 


I. INTRODUCTION 

The stochastic projected Gross-Pitaevskii equation 
(SPGPE) [1] was derived via a rigorous treatment of the 
reservoir interaction process governing the evolution of 
high-temperature Bose-Einstein condensates (BECs). The 
theory emerged as a synthesis of quantum kinetic theory [2-4] 
and the projected Gross-Pitaevskii equation (PGPE) [5-7]. 
The former gives a treatment of condensate growth in exper¬ 
iments from a particle kinetics point of view [8, 9], and the 
latter generalizes the Gross-Pitaeveskii equation to describe 
systems where the condensate is one of many degenerate 
modes — precisely the physical situation near the BEG 
transition [10-13]. The central principle allowing unification 
of the two approaches is the use of a projector to separate 
the complete system into subsystems that may be treated in 
different levels of approximation. The projector implements 
a high-energy cutoff, enabling a dynamical description of the 
low-energy system of interest, both formally [1] and numeri¬ 
cally [14-16], in the Wigner phase-space representation. This 
general procedure is familiar from the open systems theory of 
quantum optics [17, 18], and leads to a stochastic differential 
equation governing the system evolution containing two kinds 
of dissipative terms: those associated with collisions causing 
particle exchange with the reservoir (growth/loss), and those 
causing damping and diffusion of energy without particle 
transfer [19]. The simplified SPGPE derived by neglecting 
the energy damping terms (the simple growth SPGPE ) has 
been used to study the growth of BEG from a rotating thermal 
cloud [14], the formation of spontaneous vortices via the 
Kibble-Zurek mechanism [20], the decay of a vortex [15, 16], 
to the formation of a persistent current [21], to the onset of 
quasicondensation in an elongated 3D system [22], and to 
modelling thermal fluctuations in a continuously pumped 
atom laser [23]. In application to thermometry in ID, this 
description gave excellent agreement with exact Yang-Yang 
thermomety [24]. 

The energy damping collisions between low energy atoms 
and reservoir atoms are a source of phase noise [25-29], stem 
from the quantum brownian motion master equation [17], and 
cause significant damping process for highly non-equilibrium 


dynamics [19, 30]. Even in quasi-static systems, the energy 
damping process can set up superfluid counterflow if the reser¬ 
voir develops a temperature gradient [31]. Work on the full 
SPGPE that includes these terms has focused on numerically 
implementing the energy damping (scattering) terms in the 
SPGPE [19, 32], and on generalising the theory to spinor and 
multicomponent systems [33]. Despite these developments, 
there remains much to be understood about the theory and its 
implications. In particular, a complete picture of the dissi¬ 
pative dynamics of excitations in BEG such as vortices and 
solitons [15, 34, 35] requires a deeper understanding of the 
energy-damping terms. Energy damping may also play a role 
in the Kibble-Zurek mechanism of spontaneous vortex forma¬ 
tion during the BEG transition [20, 36-38], and in 2D quan¬ 
tum turbulence (2DQT) where many vortices confined to pla¬ 
nar motion interact to collectively transport energy to large 
scales [39-41]. 

Solving the three-dimensional SPGPE [19] is in general a 
numerically challenging task, due to the need for an exact 
energy cutoff (requiring a single-particle basis for numerical 
propagation [14]), and the additional expense of evaluating 
the non-local energy-damping term and integrating a stochas¬ 
tic differential equation containing multiplicative noise [32]. 
However, for systems that are tightly constrained in one or 
two dimensions, the implementation of the energy cutoff via 
a projector provides a natural route to obtaining an effective 
theory in the low-dimensional subspace of weakly confined 
dimensions. Experimentally, a high level of control over con¬ 
finement geometry is accessible optically [42-44], allowing 
the study of dissipative superfluid phenomena in new settings. 

In this work we consider systems with strong confinement 
along one or two spatial dimensions and project the 3D- 
SPGPE onto the ground state of the harmonic trap in each of 
the dimensions of strong confinement, to obtain an effective 
SPGPE in the subspace to which dynamics are confined. This 
procedure involves choosing the energy cutoff so as to elimi¬ 
nate all transverse degrees of freedom. Identifying the regime 
of validity of this procedure requires some care, as the reser¬ 
voir interaction theory we use requires that the system has a 
3D thermal reservoir; if the reservoir is also low-dimensional, 
the collision integrals determining the reservoir interaction 
rates become nontrivially spatially dependent [1, 14, 33]. We 
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assess the validity of dimensional reduction using variational 
analysis, and for the ID regime, by comparison with 3D nu¬ 
merical simulations. 


II. HIGH TEMPERATURE C-FIELD THEORY 

The projected open systems theory of high-temperature 
EEC begins by setting up an energy cutoff in an appropri¬ 
ate single-particle basis of modes representing the system. 
Where possible the basis is chosen to diagonalise all time- 
independent linear terms in the Hamiltonian, since it is then 
a good basis for imposing an energy cutoff for the interacting 
system (provided the cutoff is chosen sufficiently high). 

The system of interacting bosons confined by external trap¬ 
ping potential VrCr, t) = Vo(r) + dy(r, t) is described by the 
Hamiltonian 

d^r^^(r)['H(r) + SV(r,t)] 'P(r) 

, 1 ) 

where the time-independent single-particle Hamiltonian den¬ 
sity is 

'H(r) = -^ + Vo(r), (2) 

2m 

and u = Anfisa!m is the interaction parameter for s-wave scat¬ 
tering length a. The field operators satisfy the usual equal¬ 
time Bose commutation relations, with non-vanishing com¬ 
mutator 

mrX'P^Cr')] =<5(r-r'). (3) 

The system is represented in terms of the single-particle eigen¬ 
functions 0„(r), that satisfy 

"Hiftipnir) = e„<^„(r), (4) 

where the index n denotes all quantum numbers defining a 
unique eigenfunction. The field operator expansion then reads 

4'(r) = y]a„^„(r), (5) 

n 

for single-mode operators dn satisfying Our 

first task is to consistently separate the system into a C-region, 
where populations are appreciable and the atoms are at least 
partially coherent, and an /-region where populations are low 
and atoms are incoherent. We then seek an equation of motion 
for the C-region, treating the /-region as an incoherent reser¬ 
voir. A significant advantage of effecting the separation in the 
single-particle basis is that at sufficiently high energy it diag- 
onalises the many-body problem, thus providing a good basis 
for separating the system, provided the separation energy is 
large. We define the C-region as C = {6„ < 6cut}, where 6cut 
will be significantly larger than the system chemical potential 


p (of order 2p - 3p). In this basis we introduce the orthogonal 
projection operators 

p=Yj (6) 

Q=\-P, (7) 

satisfying PP = P, QQ = Q, Qp = 0. In the position rep¬ 
resentation the field operator decomposes into C-region and 
/-region operators as 

'F(r) = P^(r) + Q^(r) = -r ^(r) (8) 

respectively, where 

Mr) = rnr) = E Mr) f dh C(r)4'(r) (9) 

e«<ecut ^ 

= 2] ^nMr), (10) 

defines the spatial representation of P and the projected field 
operator with commutator 

[.A(r), ^^(r')] = <5(r, r') = Ur)^^). d D 

This formal separation of the system provides a natural ap¬ 
proach to deriving an equation of motion describing the evo¬ 
lution of the C-region, the details of which can be found else¬ 
where [1, 33]. The derivation proceeds by mapping the mas¬ 
ter equation for the C-region density operator to an equation 
for the Wigner distribution of the system; an equivalent diffu¬ 
sion process for a classical field j/r(r) (c-field) is then obtained, 
the moments of which correspond to symmetrically ordered 
averages of the field operator at equal times [30, 45]. The 
/-region field pir) = (3'F(r) enters the equation through the 
rate functions determining the strength of reservoir interaction 
processes [14, 19]. 

For our purposes, we take as our starting point the three- 
dimensional Stratonovich SPGPE for the C-field [19]. Taking 
our energy reference as p, the SPGPE takes the form 

(5')Jj/^(r, t) = , (12) 

with 

ihdil/\ =P{Lil/dt}, (13) 

\H 

= P {-iyLj/dt ifidWir, t )}, (14) 

{S)indil/\ = P {y(r, t)j/dt - nif/dUir, t )}. (15) 

\e 

The Hamiltonian evolution relative to the reservoir chemical 
potential is generated by 

Lt//(r, 1) = ['H(r) + 6Vir, t) + Op - yu] Mr, t), (16) 

thus recovering the PGPE [7] 
dil/ir, t) 

( 17 ) 

ot 
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This equation of motion describes the evolution of a low- 
energy fraction of atoms with partial coherence, including the 
condensate and many excitations [7, 46, 47] [48]. The high en¬ 
ergy reservoir coupled to the PGPE is described by chemical 
potential (//), temperature (T), and cutoff energy (6cut), and the 
functions 


8^2 ~ 

^ “ J2 2-i e2/3ecu.y ® 

dB j=l 


. 1.7 


^2/5 6c 


Jd^r's(r-r')V' 


y(r, t) = -n 

Ifl 

j(r,0 = ^[W-rv^^]> 


e(r) = 
A\ = 


r , 


M 
(2;r)3 

16710^ 

^/5(^CUt~/^) - t 


(18) 

(19) 

( 20 ) 
( 21 ) 
( 22 ) 


where AdB = , (3 = l/^sT, and 0[z, x, a] = 

2^0 ^^/(^ + Is the Lerch transcendent, and where y and 
s(r) are both dimensionless. The noise terms are Gaussian, 
with non-vanishing correlations 


{dW\r, t)dW(r\t)) 
{dU(r, t)dU(r',t)} 


n 

(23) 

n 

(24) 


For a reservoir close to equilibrium, the growth rate G(r) 
(See Eq. (6) of Ref. [19]) acquires the same value over most of 
the C-field region, and varies slowly near its edges [14], justi¬ 
fying our approximation of the growth in the 3D system by its 
spatially independent form, and defining the growth parameter 
fiG(r)/kBT ^ y. This approximation has the additional advan¬ 
tage of putting the treatment of growth and scattering on an 
equal footing: for a thermal cloud near equilibrium the coef¬ 
ficient of the scattering term is spatially invariant (aside from 
the irreducibly non-local scattering kernel). In (12) and what 
follows we use the subscript s (y) to denote the the scattering 
(growth) reservoir interaction terms in the time-independent 
reservoir approximation used in this work. We emphasize that 
the energy cutoff defining the reservoir is taken at quite a high 
energy (relative to ja), and the dimensionless damping rates 
are quite small (~ 10“"^ - 10“^ in experiments [21]), so that in 
practice this approximation allows for significant dynamics of 
the non-condensed atoms in the C-field. 

Defining, respectively, the C-field atom number, Hamilto¬ 
nian, and grand-canonical energy as 


N = 


(25) 

H = 


(26) 

K = 

H-pN, 

(27) 


we can get some insight as to the role of the terms in the equa¬ 
tion of motion. Without noise, i.e. formally setting all noises 



FIG. 1. (Color online) Schematic of Low-D SPGPE regimes for 
an anisotropic 3D Bose gas. (a) Reservoir interaction processes for 
number damping (y) and energy damping (e) respectively, (b) ID 
and 2D scenarios of anisotropy with low-D c-field and 3D thermal 
cloud, (c) Separation of energy scales due to anisotropy allows the 
transverse degrees of freedom to be treated as part of the thermal 
reservoir. 


to zero in the SPGPE, the equations of motion generate the 
dissipative evolution 

dN 2y 

— = (28) 

dt n 

(29) 

where ii{t) = f d?Y il/*(L -r ii)if/IN is the instantaneous chemi¬ 
cal potential (the definition (16) introduces a shift by yu), and 

f = -||^=rWP-»A<|:^|k.j(k)P (30) 

is non-negative. Clearly y describes number damping (with 
associated energy damping), while s describes energy damp¬ 
ing (without associated number damping), and both processes 
arise from scattering between C-region and /-region atoms, 
with rates that are quadratic in the scattering length. The two 
collision processes responsible for stochastic evolution of the 
SPGPE are shown schematically in Fig. 1 (a). 


III. DIMENSIONAL REDUCTION 

To introduce notation, we first consider the reduction to a 
2D subspace, i.e. a system with tight parabolic confinement 
along the z-axis, defined by harmonic oscillator frequency 
Defining = (-x,}^), we then have 


moj^z 

yo(r)^yo(rJ+^^. (31) 

We can ensure that the system has a 3D thermal cloud by im¬ 
posing the temperature condition 

HclJz < ksT. (32) 

Then it will be a good approximation to introduce the en¬ 
ergy cutoff as 6cut ~ fioJz/2, i.e. we treat modes above the 
ground state in the z-dimension as part of the /-region. Since 
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the energy scale in the z-direction is quite large, the mode 
number in the 2D subspace can still be varied, as 6cut can be 
chosen in the range ^6 g)^/ 2 < 6cut < 3lz6c)^/2 while still main¬ 
taining a formal dimensional reduction. The dynamics of the 
superfluid will remain confined to the plane provided the 2D 
chemical potential satisfies jj 2 ^ The validity condition 
for dimensional reduction can be expressed as 

jd2 ^ Jiojz ^ ^bT. (33) 


The situation is very similar in ID. The trapping potential 
is now assumed to have strong parabolic confinement along 
the y- and z-axes: 


VoCr) ^ yo(v) + 


-r z^) 

2 


(34) 


Putting the two conditions together, a consistent 3D reservoir 
theory of the ID subsystem is obtained provided the condition 


yUl < ksT 


(35) 


is satisfied. It will be helpful in what follows to define the 
D-dimensional transverse frequency as 


Qd - 



(D = 2), 
(D = 1), 


in terms of which our validity condition reads 


(36) 


fij) fi^j) < kBT. 


(37) 


These two scenarios of anisotropy are shown schematically 
in Fig. 1 (b); the separation of energy scales is illustrated for 
the pancake configuration in Fig. 1 (c). The first inequality 
in (37) raises a central question of this paper: how different 
must these scales be for validity to be ensured? We will return 
to this question in Secs. IV and V. In what follows we will 
also find it convenient to introduce the oscillator length scale 
associated with the transverse trapping potential: 


To keep notation simple, we will often refer simply to Q and 
(T as it will be clear from the context which quantity is being 
used. For parabolic traps it also customary to define the oscil¬ 
lator lengths ax = ^Jfilmojx, etc, and these quantities will be 
reverted to where appropriate. 

In each of the tightly confined dimensions we insist that the 
C-field is separable, and that the system is in the transverse 
ground state wave function, e.g. 

/ 1 

0o(z) = (39) 

where the subscript of ctd will be omitted unless it is essential 
for clarity. Dimensional reduction is carried out by integrating 
over the tightly confined dimensions. We can state this pro¬ 
cedure formally noting that we can always (since our basis is 
separable) write the projector in the factorized form 


P = PdP3-d, 


(40) 


where is the reduced dimensional projector for the D- 
dimensional subspace [the low-D form of (9)], and P^-d is 
the projection onto the transverse degrees of freedom; for the 
regime we consider, the transverse projector can be evaluated 
explicitly as the projection onto (39) in each dimension. 

To avoid overly cumbersome notation we denote the dimen¬ 
sionality of the wavefunction by its spatial argument r, r^, and 
X in 3D, 2D and ID respectively. 


A. Number damping 


Since the number damping rate y is well approximated as 
spatially invariant, the number damping terms in low-D are 
formally identical to their 3D form. Dimensional reduction 
leads to a modified interaction strength, as has been described 
in many previous works on the GPE [45]. The dimensionally 
reduced forms for V, H, K may be summarized as 


^D = f d°v \^|,\\ (41) 

Hd = ^d’^v r (^D + y (42) 

KD = HD-^lDND, (43) 

where the interaction parameter 

{Anfisa! m (D = 3), 

gz) = I aimer (D = 2), (44) 

[inna, (D = l\ 


is obtained by integrating out the transverse ground-state 
wavefunction in 2D and ID respectively, and the chemical po¬ 
tential is given by 


I^D 


// (D = 3), 

^ Id - fi(jJzl2 (D = 2), 

fi - (D = 1). 


The single-particle Hamiltonian density is 




"Hd - 


< -n^V\l2m + VT{Y^j) 
-fi^dl/2m Vt(x, t) 


(D = 3), 
(D = 2), 
(D = 1). 


(45) 


(46) 


The dimensional reduction applied to (13) and (14) then gives 


ifidif/ = Pd {(1 - iy)LDil/dt iMWoiY, t )}, (47) 


where 


Ld^ = (TTz) -I- goW? - (48) 


(JW* (r, t)dWD(r\ t)) = v)dt, (49) 

n 

dz)(r',r) is the D-dimensional projected delta function, and 
r takes the form appropriate for subspace D. When conve¬ 
nient, we will also refer to these projectors as Px and P^ in 
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ID and 2D respectively. Thus, in the regime where the ther¬ 
mal cloud remains 3D (i.e. well away from the extreme low-D 
limit where S-wave scattering is modified), the dimensional 
reduction for the number-damping terms yields only a renor¬ 
malized interaction strength, and a new reference energy for 
the chemical potential. 


B. Energy damping in ID 


We will describe one dimension derivation in detail to give 
an indication of the procedure for formally obtaining the equa¬ 
tion of motion for the theory in a reduced dimensional sub¬ 
space. The energy damping terms are found by projecting the 
full field j/r(r, t) = 0o(}^)0ofe)<A(v, t) as 


(S)ifidi//(x, o| dy dz ^o(y)<Po(z)Pj, 
y(r, t)dt 

- fnl/{Y, t)dU(Y, ok 


This gives the one-dimensional form 

{S )ifidi//{x, o| = (v, t)if/(x, t)dt 

- fiil/(x, t)dUi(x, t) 


(50) 


(51) 


with new potential and noise obtained by evaluating the pro¬ 
jectors for each term. We thus define 


Vi{x,t)= J'^yJ' dz\(po(y)\^\(pQ{z)\^V{r,f) 

~ J' dx' s\(x - x')dx' j(x', t). 


(52) 


in terms of the function 

M 


Si(x - x) = 


(27r)3 


J dkx J 


d^kj 

1 ^ 




d2 g<k,.(r,-rl)! 


-rl/cd 


M 


(2n')^cr 


f r 

d do Ju- 


du 


^Ju^ (kxCr)^ 


(53) 


where we have changed variables to u = kj_cr. Evaluating the 
integral, we obtain 


where 


^i(-^) = t;— I dke^^^Si(k), 

271 J-oo 




I V2/ 


(54) 


(55) 


2 

and G(q) = erfc(^) is the scaled complementary error func¬ 

tion. 

Projecting the noise onto the ID subspace, we have 


dU\{x, t) = 



dz \(l)Q{y)^\(l)Q{z)\^dU{r,f). 


(56) 


Since the noises are completely determined by their first and 
second moments, we can instead evaluate the correlation func¬ 
tion 


{dUfx, t)dUi{x\t)) 


dz \Uy)nUz)V 


X f dy'J' Jz'|^o(>’')Pl0o(z')P 
x(dU(r,t)dU(r',t)}, (57) 


and evaluating the integrals gives an expression for the noise 
correlator in terms of (54) as 

{dU\{xj)dU\{x J)) = — ^£i{x - x)dt. (58) 
n 

This reduced correlation function gives a complete description 
of the noise for the ID subspace. 


C. Energy damping in 2D 


The energy damping terms are projected from the full field 
</6(r, t) = j/r(r^, t)(f)o(z) as 

(S)iMiJf(r^,t)l= J dzmPxy\^ 

y(r, t)il/{Y, t)dt 

-rulj{r,t)dU{r,t^ (59) 
giving the two-dimensional form 

(S)i%dil/(r_i_, f)| = Pi\ 0<A(r-L> t)dt 


nipir^, t)dU2{Y^, t) 


(60) 


by evaluating the projectors on each term. We thus define 


V2(r 




dz \4>o(z)\ V(r,t) 


-nj dVl£2(ri-rl)Vl.j^(rl,0, 


(61) 


where 

e2ir± - rl) ^ 


M 

(27r)3 


ik-iz-z') £_ 


r f dk,-l 

*J d—oo 

fdzf dz' 

r°° d 

J-oo -Ju 


oirx-r'x) 

du 




(62) 
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Evaluating the integral, we arrive at the 2D result 

£2(r^ - rl,0 = ^ J (63) 


where 




(64) 


and F(x) = e^K^ix) is the scaled modified Bessel function. 
The noise is 


dU2{r_ 




dz \<Pq{z)\ dU{Y,l), 


(65) 


and as above, we evaluate the noise correlation in 2D 

{dU2{r^,t)dU2{r'^,t)) = J dzIMzff dz'IMz'f 

x(dU(r,t)dU(r',t)), (66) 


SD(r) = ^ f d°k e'^^z)(k), (72) 

and noise correlation functions 

(dW^(r, OdWoir', t)) = ^r<iz)(r', r)Jr, (73) 

n 

{dUoix, t)dUD{Y\ t)) = ed(x - r')dt, (74) 
n 

where 

= (75) 

n 

is the projected delta-function of the D-dimensional C-field 
region. We have used the shorthand 0i^\r) to represent the 
modes in the D-dimensional subspace, i.e. for a harmonic 
oscillator basis in D = 2, the modes are a separable prod¬ 
uct of independent harmonic oscillator eigenstates: 0i^\r) = 
fnM)^nyiy)> 


finding 


{dU2(r^j)dU2(r'^J)) 


2kBT 

n 


£ 2 (r^ - r'Jdt. 


(67) 


D. Summary of Low-D SPGPE 


In this subsection we summarize the SPGPE in each of the 
three regimes. To simplify notation hereafter we drop the 
cumbersome subscript r^, as it will be clear form the con¬ 
text that the wave function, spatial coordinate, and wave vec¬ 
tor coordinates are set by the dimension D. We define the 
Eourier-space scattering kernel 




ik|-i 


(D = 3), 
(D = 2), 
(D = 1), 


( 68 ) 


The asymptotic scaling of S z)(k) is 




(3-D)/2 


(69) 


E. Separation of length scales in ID: White noise limit 


While the expressions (54), (55), and (58) give a formally 
exact reduction to ID, a further approximation can be identi¬ 
fied in the regime of high anisotropy where the length scale 
cr becomes much smaller than the length scales over which 
</r(x, t) changes. The scattering kernel may then be approxi¬ 
mated as 

5i(fc)^5i(0)= —(76) 

Vs^ 


giving the potential 


Vi (x, t) ^ Vi (v, t) = -fiMd^fixf (77) 


with 


M = 


M 

VStt^ 


(78) 


Making use of the separation of scales to simplify the noise, 
we arrive at 


£i(x) ^ £i(x) = Md{x), (79) 


where for D < 3 there is enhancement by due to con¬ 

finement. 

The low-D SPGPE may be summarised in the form 

{S)ifidil/{Y, t) = !Pd|(1 - t)dt + iftdWuix^ 0 

-r Voix, 0<A(r, t)dt - fiil/ix, t)dUD{^, o|, 

(70) 

with potential 

Voir, t) = -n J d°r'so(r - r')V' ■ jz,(r'), ( 71 ) 


and the white noise correlation function 

(dUi(x, t)dUi(x, 0) ~ (dUi(x, t)dUi(x, t)) 

= M6{x - x')dt. (80) 
n 

Thus, for large enough anisotropy the energy-damping noise 
becomes broadband white noise. 


IV. VARIATIONAL ANALYSIS 

In this section we present an assessment of the validity cri¬ 
teria (33), (35) for the effective low-dimensional stochastic 
equations of motion. 
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A. Estimating SPGPE parameters 


A basic problem in SPGPE theory involves determining pa¬ 
rameters for modelling experiments, as has been studied in 
3D [15]. In general we assume the experimental parameters 
take the form of a definite total particle number N, and tem¬ 
perature r, and we wish to find /j(N, T) and T) suitable 

for modelling the system using SPGPE theory. Eor the low-D 
effective theories, the value of 6cut is partially constrained by 
the geometry, requiring that we assess the validity of project¬ 
ing into the low-D subspace. 

To determine system parameters, we start from an ideal gas 
description of the thermal cloud in a 3D parabolic trap, with 
BEG transition occurring at temperature 




N 

1^)1 


(81) 


for a given total particle number N, and geometric-mean trap 
frequency = oJxOJyOJz. Eor highly anisotropic confinement 
we must also include the finite-size and interaction shifts in 
the determination of giving the transition temperature 


B. Corrections to the transverse ground state: hybrid 
Lagrangian variational method 

The low-D chemical potentials given by (86) describe the 
regime where the transverse wavefunction is in the ground 
state. The regime of validity for this approximation can be 
understood by determining the leading order correction to this 
expression, for a transverse wavefunction that is allowed to 
vary, while remaining Gaussian. A variational treatment of 
dynamics with such a transverse wavefunction leads to the 
hybrid Lagrangian variational method (HLVM) [50], giving 
an effective Gross-Pitaevskii equation with modified interac¬ 
tions, coupled to an auxiliary equation of motion for the trans¬ 
verse width Aj). This description provides an estimate of the 
importance of post ground-state corrections in ID and 2D. 

The HLVM method seeks to find a low dimensional effec¬ 
tive GPE, starting from the 3D GPE Lagrangian density 

X = ^ ((D(D; - (D,(D*) + ^|V(D|2 + yKr)|<Dp + ^|<D|^ 

2 2m 2 

(89) 


Tc(N) = T° + STf + dry. 

(82) 

Applying the variational principle to (89), 

one obtains the 

for finite-size and interactions shifts [49] 


Euler-Lagrange equation 


cy^fs _ ^,^-1/3^0 

2^(3)2/3a> 

(83) 

V a = 

Z-i ''Is®*/ <90*’ 

r]=x,y,z,t ' ” ' 

(90) 

dry = -1.33riV‘^®T°. 

(84) 



Here oj = (ojx ojy (0^)13, and a"^ = a^OyO^ is the geometric 
mean of the oscillator lengths Ox = ^|hJm^o~x, etc. Eor a given 
atom number A, we estimate the condensate number via the 
ideal gas equation of state 


Ao/A = 1 - (r/r,)l 


(85) 


Eor Ao atoms in a low-D Thomas-Eermi state, we use (44) 
to find the chemical potentials 




(^ = 2), 

We then use (45) to estimate the 3D chemical potential as 


( 86 ) 


K + (D=l). 


(87) 


which leads to the 3D GPE as equation of motion for the field 
O. The basic approach in Ref. [50], extended here io D = 1, 
is to assume a trial wavefunction of the form 


O(r,0 


_ |</^(v, y, t)A{t) exp [-z^ j2A^ + iKz^^ (D = 2), 
ykiz, t)A(t) exp i^-r\l2A? + {D = 1), 


(91) 


describing the lowest energy (breathing) excitation in the 
transverse dimensions. Integrating out the transverse dimen¬ 
sions, we then find the Lagrangian densities 

-Cd = -Jr + kz)(r, 

2 2m 2 


+ 0d(OI<AI 5 


where gives the Laplacian in D dimensions. 


(92) 


We determine the D-dimensional energy cutoff by inverting 
the Bose-Einstein distribution 


= log 


1 




IcbT /2d, 


( 88 ) 


goit) - 


(27iA(t)^)(^-^y^' 
is the effective interaction, and 


(93) 


where the population of modes at the cutoff energy and typi¬ 
cally chosen to be n^^^ ~ 1 - 3, consistent with the truncated 
Wigner approximation [30]. This choice of 6^^ ensures that 
the number of modes in the weakly confined direction is pre¬ 
served when the dimensional reduction is made. 


0Z)(O 



4mA^ m 4 

, in^A^K^ , mcojA^ 
2mA^ m 2 


(D = 2), 
(D = 1). 


(94) 


contains all of the transverse energy terms. The Euler- 
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Lagrange equations are now 


dt 


dZo] 

1 _ 

dk j 

4 Ok 

d£D\ 


dX 1 

dA 

dZ ) 

1 dZ 


~ dr' 


(95) 

(96) 

(97) 


where 77 ' runs over t and the spatial coordinates of the D- 
dimensional subspace. For either dimension, using (95), inte¬ 
grating the result J^r, and using dtN = 0, gives k = mA/2fiA; 
this expression allows the simplification of the A equation 
from (96), as in Ref. [50]. The end result of this procedure 
is the set of coupled equations of motion 



A — 


gpiPCpit) 

rri^A^ AmN 


(98) 

(99) 


where 


Coit) = 



( 100 ) 


and to arrive at Eq. (98) we we have transformed the wave- 
function to the rotating frame i}/ = if/ exp (-/ The 

time-independent equations 


quadratic equation - fi^ I= 0 , with so¬ 

lution 


4 = 4 


1 

12 ^ 0 ) 


1/2 


+ CTn 


aDli% 
^ inkin' 


(106) 


where Qo and ctd are defined in (36), and (38). Thus, to lead¬ 
ing order 







(107) 


and the low-dimensional regime will be reached when the sec¬ 
ond term is much less than unity. We can thus limit the frac¬ 
tional change in cr^ to be less than S by imposing the restric¬ 
tion 


Md ^ 


2fLQ.j) 

- o. 

ao 


(108) 


For parabolic trapping, a choice of d = 0.2 translates to the 
constraints yUi < fioj^l2 and 7/2 < fioJzl^, and as we shall see 
in the next section, the boundary of the ID regime is well 
approximated by this choice. In general the specific trapping 
geometry of the low-D space must be considered to determine 
ao- 


V. ID REGIME OF VALIDITY 




2m 


+ 14)(r) + 


qIa = 


_i_ Sd^d 
m^A^ AmN ’ 


( 101 ) 

( 102 ) 


can now be solved in the Thomas-Fermi approximation for 
N = No atoms in a Thomas-Fermi condensate with chemical 
potential For a parabolic trap, we make use of ( 86 ) (with 
(To replaced with variational parameter Au) to obtain expres¬ 
sions for A/q, and find the Thomas-Fermi expressions 


{ {Inf 
_ I moj\ 
“ I 16 


(D = 2), 
(D = 1). 


We can then obtain the simple result 


gpCp 

No 


aD^^D■ 


(103) 


(104) 


For the case of parabolic confinement, the trap and geometry- 
dependent numerical factor is 


In this section we test the ID validity criteria by considering 
the equilibrium state of a finite temperature Bose gas in a cigar 
trap, as found from numerical simulation of the ID-SPGPE, 
and compare the results with the HLVM criteria (108), and 
with the 3D-SPGPE. 


A. ID System 

In order to test the validity of our low-D description, we 
consider the dimensional reduction of a cigar-shaped trap. 
The trapping potential of the system is given by 

^(r) = y + 2^)) > (109) 

where ojx ^ 0J±. We use a system consisting of 1 x 10^ 
total number of ^^Rb atoms held in a trap with constant ojx = 
2n Hz. In terms of the C-region and /-region, this is comprised 
oiN - Np Ni. To investigate the transition towards the ID 
regime, we vary via the aspect ratio parameter 


|2/3 (D = 2), 
(4/5 (D=l). 


(105) 


. _ Clz 

2 V = — —, 

aj_ cot 


( 110 ) 


Using (104) in (102), and denoting the stationary solution so that as A increases the system becomes increasingly pro- 
for Ap in Thomas-Fermi approximation by Ap, we find the late. 
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FIG. 2. Phase plots demonstrating the regime of validity for the one¬ 
dimensional SPGPE theory, for a range of relative temperatures and 
length aspect ratios (A). Both the total atom number (N = 1 x 10"^) 
and ojx = 27 t Hz are fixed for all T and A, and jui is calculated via 
( 86 ). 


B. SPGPE regimes 


In this section we highlight the different regimes which oc¬ 
cur within the SPGPE theory during the transition to the ID 
regime. We calculate the chemical potential via Eq. (86), 
for 1 < A < 50, 0.5 < TjTc < 1. Different regimes 
of validity arise from considering the ID validity criteria in¬ 
equalities in Eq. (35). The key energy ratios are shown in 
Eig. 2. We firstly consider and see that this ratio de¬ 

creases with increasing anisotropy, and most dramatically at 
higher temperatures. In contrast, Ttco^/kET increases with in¬ 
creased anisotropy, eventually suppressing the transverse ther¬ 
mal cloud, for large A. In order to make the dimensional 
reduction within the SPGPE framework, the system must be 
anisotropic enough so that the transverse wavefunction is in 
the ground state, but still in a regime where the thermal cloud 
is three-dimensional so that the SPGPE is valid. 

An important consideration is the role of phase fluctuations 
when approaching the quasi-ID regime. The characteristic 


A 

4 

8.5 

14 

24 

idilho)^ 

0.83 

0.50 

0.36 

0.25 

ftcxJ^^lksT 

0.18 

0.31 

0.45 

0.70 

^cut / P 

2.25 

1.88 

1.66 

1.47 

«cut [3D] 

1.87 

1.99 

1.80 

1.82 

^CUt [ID] 

1.90 

1.97 

2.04 

1.99 

Ni X 10-3 

7.57 

5.93 

5.90 

6.05 

Ni X 10-3 

5.80 

5.86 

5.90 

6.06 

No X 10-3 [3J3] 

6.22 

2.13 

0.75 

0.22 

No X 10-3 [ID] 

4.7 

2.35 

0.83 

0.23 

TIT^ 

0.12 

0.85 

3.07 

11.95 


TABLE I. Summary of system parameters used for ID and 3D 
SPGPE simulations. Equilibrium properties are also shown. For 
all values of A we have constant = 2?: Hz, A = 1 x 10"^, and 
TITc = 0.75. 


temperature for the onset of phase fluctuations is [51] 

T^ = N{noj,flnkB, ( 111 ) 

which is shown in Fig. 2 via the ratio the system chem¬ 
ical potential yu, is approximated as yUj defined in (86). If the 
system is in the phase fluctuating regime, the condensate pop¬ 
ulation is suppressed below the 3D ideal or interacting gas 
results. This does not alter the validity of the ID-SPGPE 
regime, but is important to consider when interpreting the re¬ 
sults of simulations, or in formulating approximate treatments 
(such as variational approaches) that require a coherent wave- 
function. 


C. Transition to the ID regime - SPGPE simulations 

We now quantitatively investigate the transition from the 
3D to the ID regime. We calculate equilibrium states of the 
SPGPE for fixed TjTc = 0.75, and 4 < A < 24. Full 
simulation parameters are shown in Table I. We And equi¬ 
librium states of the SPGPE by evolving the simple-growth 
SPGPE with y = 0.4, so that finding equilibrium only re¬ 
quires that we evolve the system for t = A trap cycles of evo¬ 
lution. To compare between the one-dimensional and three- 
dimensional equilibrium states, we consider the ensemble av¬ 
eraged classical-field density ^(r) = (|0^(r)p), and the conden¬ 
sate density ;^(r) found via the Penrose-Onsager criterion: 

J dr'p(r, v')x(r') = Nox(r), ( 112 ) 

where the one body density matrix p(r, r') = (ik(r)ik'^(r')) is 
constructed from an ensemble of 2000 trajectories for each 
parameter set, and the largest eigenvalue No gives the con¬ 
densate occupation. Note that the ideal gas estimate for the 
condensate population is Aq = 5.7 x 10^ for all values of A, 
and our system parameters are chosen so that finite-size and 
interaction shifts to Tc are taken into account [see (82)], to 
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x/xo 



FIG. 3. (Color online) Ensemble averages [using 2000 trajectories] 
of the classical-field density profile nnix), for a range of anisotropies 
at r = 0.75Tc, comparing the effective one-dimensional implemen¬ 
tation (red) with three-dimensional simulations (blue). The density 
corresponding to the condensate Xd{^) is shown by the dashed red 
and blue lines. The ID Thomas-Fermi density for chemical potential 
jui is shown for comparison (black dash). 


ensure that TjTc stays constant (as does + A^/) as A 

increases. 

To quantitatively compare three-dimensional and one¬ 
dimensional densities, we integrate out the tightly confined 
dimensions of the three-dimensional density, i.e. 


nsW = 



dz{W{x,y,z)\^). 


(113) 


This projection extracts the correct ID density from the 3D 
theory, for comparison with our effective low-D theory. 

In Fig. 3 we compare the total c-field and condensate den¬ 
sities, for the 3D and ID simulations. For the smallest aspect 
ratio, A = 4, we find that there is a significant difference be¬ 
tween the 3D and ID density profiles. The average c-field 



FIG. 4. (Color online) Summary of the different regimes for a quasi- 
ID system as a function of TjTc and A. The regimes as described in 
the main text are separated by the lines correspond to = 1 (red 
line), fico^lksT - 1 (blue line), and jiilficx)^ = 1/2 (thick grey line). 
The circles correspond to SPGPE parameters used in the simulations 
of section V C [See Table I]. 


number Nc, and condensate number Aq, are shown in Table I, 
and for this aspect ratio both the ID-SPGPE values are ~ 24% 
smaller than those of the 3D-SPGPE. The large discrepancy 
stems from the inconsistency of the dimensional reduction: 
the choice of cutoff energy includes some transverse modes in 
the c-field, thus incorporating modes of the reservoir into the 
ID-SPGPE description. Increasing the aspect ratio to A = 8.5 
— the value that ensures that the inequality (108) holds — the 
population artefact vanishes, and we find excellent agreement 
between particle densities in the ID- and 3D-SPGPE. The 
condensate mode density is also closely comparable, and the 
c-field density agrees quite well with the ID Thomas-Eermi 
result. The ratio = 0.85 < 1, and the condensate popula¬ 
tion is suppressed below the ideal gas estimate, but the ID and 
3D-SPGPE values for Nq are in close agreement. Eor A = 14 
the condensate population shows an order of magnitude depar¬ 
ture from the 3D ideal gas result expected for a phase coherent 
system, consistent with = 3.07 » 1. The ID-SPGPE 
describes this situation reliably, as seen in Eig. 3 where the 3D 
and ID condensate densities agree fairly well given the finite 
ensemble. Eor A = 24 the system is in the phase fluctuating 
regime ~ 12, the condensate is strongly suppressed, and 
the thermal tails of the total density are a dominant feature 
of the distirbution. Note that for all values of A > 8.5 the 
3D- and ID-SPGPE give quite similar results for Nc and Aq, 
as shown in Table I. We have also checked the cutoff depen¬ 
dence of the ID-SPGPE, as described in the Appendix, where 
we find that the cutoff can be chosen consistently within the 
classical field approximation [30] by setting the cutoff mode 
population ^cut ~ 2. 
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FIG. 5. (Color online) Average c-field density of ID-SPGPE equi¬ 
librium states, comparing different mode populations at the cutoff 
energy, for A = 8.5 and 14. 


D. Summary and Discussion 

The validity criteria (32) and (108), and the results of our 
simulations are summarized in Fig. 4, where we may identify 
several regimes: 

I. 3D SPGPE. — The system contains a 3D thermal cloud 
since %oj^ < ksT. The superfluid is also in a 3D regime with 

> 1/2, and a well defined condensate exists since 
T < T(f,. Thus this regime requires the 3D SPGPE for its 
description. 

II. Phase coherent ID SPGPE. — Here jiilhco^ <1/2 and 
thus the superfluid is ID, and iiilhco^ < 1/2, justifying the 
dimensional reduction. However the system remains outside 
the phase fluctuating regime since T < T(p. The precise value 
of iiilhoj^ for this transition is unclear, and we represent this 
uncertainty using a thick line in Fig. 4. 

III. Phase-fluctuating ID SPGPE. — In this regime T > T(j), 
and significant phase fluctuations suppress the condensate 
population relative to what may be expected from the 3D ideal 
gas equation of state (85). The reduction to a ID SPGPE re¬ 
mains valid, as hoj^ < ksT so that the thermal cloud is 3D. 

IV. Invalid ID SPGPE. — Here we have ksT < hco^, and 
the thermal cloud loses its 3D character. The dimensionally 
reduced description derived in the present work is no longer 
valid in this regime; an SPGPE theory could be obtained by 
accounting for reservoir interactions between c-field atoms 
and atoms of a ID reservoir. We emphasize that the ID- 
SPGPE still provides a valid means of creating equilibrium 
ensembles, even though it no longer describes the real-time 
evolution of the ID system. 

V. Phase fluctuating 3D SPGPE. — The condensate is sup¬ 
pressed by phase fluctuations in this regime due to the high 
anisotropy and low temperature, despite the system satisfying 
the 3D thermal cloud criterion hoj^ < ksT, and containing a 
3D superfluid. 


Fig. 4 shows that we have a limited parameter space where 
the dimensional reduction is valid. As we are always inter¬ 
ested in a reduced dimensional c-field regime, the condition 
(108) is inflexible. However, the 3D thermal cloud condi¬ 
tion hklj) < ksT can be relaxed for certain scenarios where 
an SPGPE theory still provides a valid description of a low- 
D subsystem. The most obvious example is provided by a 
system consisting of two distinct atomic species, for which 
an SPGPE theory has been derived [33]. For a regime of 
confinement (perhaps achieved through combined magneto¬ 
optical trapping) where one component is below and the 
other is above the number damping terms are completely 
suppressed, and the SPGPE takes a form where only the en¬ 
ergy damping reservoir interaction terms occur. The degen¬ 
erate component evolves according to an SPGPE describing 
energy-exchange with the non-degenerate component, and the 
3D thermal cloud constraint may be lifted. A related sys¬ 
tem involving buffer gas interactions as a mechanism for de¬ 
coherence has been explored [52]; the master equation de¬ 
rived bears some resemblance to the Quantum Brownian Mo¬ 
tion master equation governing energy damping in the two- 
component SPGPE theory [33]. 

Finally, we discuss the relationship of the SPGPE to other 
theoretical treatments of finite-temperature Bose-Einstein 
condensate (BEG) dynamics. These fall broadly into two cat¬ 
egories: generalized mean-field theories that treat the con¬ 
densate within a symmetry breaking approach [53], and theo¬ 
ries stemming from a phase-space representation of the field 
theory [1, 54]; for reviews see [30, 55]. A notable excep¬ 
tion is given by the number conserving approach to the mean 
field theory, that avoids breaking U(l) symmetry, while de¬ 
scribing the dynamics of the condensate and its Bogoliubov 
fluctuations [56-58]. A stochastic Gross-Pitaevskii treatment 
of the high temperature BEG derived via the Keldysh ap¬ 
proach to non-equlibrium dynamics [54] has also been used 
in several studies [34, 37, 59-65]. The reservoir interaction 
in this approach arises in a finite-temperature treatment of 
the many-body T -matrix, giving a theory of the low-energy 
field where fluctuations arise through the exchange of parti¬ 
cles with the thermal fraction. This approach is essentially 
equivalent to any non-projected stochastic Gross-Pitaevskii 
equation for which the equilibrium ensemble is grand canon¬ 
ical [36, 66]. It typically gives a very good description of 
equilibrium properties, and in many situations describes the 
same physics as the simple growth SPGPE. However, the the¬ 
ory lacks a formal separation of the system via an orthogonal 
projector. Consequently, there are challenges in handling the 
inherent UV-divergence in a physically and numerically con¬ 
sistent manner (with sufficient care this can be achieved in 
equilibrium [64]), and the dynamical predictions of the theory 
differ as the energy-damping terms are absent. 


VI. CONCLUSIONS 

The 3D-SPGPE has had some success in describing exper¬ 
iments with either minimal [20] or no use of fitted parame¬ 
ters [21]. However, the 3D theory remains numerically chal- 
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FIG. 6. (Color online) Average c-field density of 3D-SPGPE equi¬ 
librium states, comparing different mode populations at the cutoff 
energy, for A = 8.5 and 14. The ID-SPGPE equilibrium solution 
is shown for comparison (red line), and in close agreement with the 
3D-SPGPE for 4^, = 2 (blue line). 


lenging, and is also an inconvenient framework for developing 
analytical models of the dissipative dynamics of excitations in 
finite-temperature BECs. In this work we have presented a 
dimensionally reduced stochastic projected Gross-Pitaevskii 
equation that provides a more tractable approach for scenar¬ 
ios where a quantum-degenerate fraction of the atoms resides 
in a physical subspace of reduced spatial dimension, and the 
tightly confined transverse dimensions are well described by 
the quantum mechanical ground state of the system (assumed 
to be Gaussian). In this regime the SPGPE projection for¬ 
malism for separating the system into coherent (C-field) and 
incoherent (/-field) regions provides a physically natural and 
formally rigorous way to find the low-D SPGPE. 

We have gone some way towards validating this approach 
for the ID-SPGPE, via a variational treatment of the trans¬ 
verse degrees of freedom and numerical simulations compar¬ 
ing 3D-SPGPE with ID-SPGPE. We find that the obvious in¬ 
equality governing the regime of validity, namely that the sys¬ 


tem chemical potential, yu, and transverse oscillator energy, 
hojj_, should satisfy yu fioj^ ksT, is overly restrictive and 
can be significantly relaxed. In practice the condition for the 
ID-SPGPE can be stated as Kfi < < ksT, where a: > 1, 

and the precise value is system dependent; for the cigar trap, 
we find empirically that k = 2 provides a reasonable guide for 
validity of dimensional reduction [see (108)]. 

Areas of future interest include quantifying the 2D-SPGPE 
validity regime, application of the low-D SPGPE to the 
Kibble-Zurek mechanism, to the dissipative dynamics of soli- 
tons and vortices, and to two-dimensional quantum turbu¬ 
lence. 
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Appendix: Cutoff dependence 

We have carried out two separate tests of the cutoff depen¬ 
dence of our results. Eirst, we varied the specific choice of 
population at the cutoff energy, for the ID-SPGPE. The re¬ 
sults are shown in Eig. 5, where we see essentially negligi¬ 
ble change in the particle density, despite changing the cutoff 
population by 33%. Second, we performed 3D-SPGPE sim¬ 
ulations for different choices of cutoff population and com¬ 
pared with the ID-SPGPE results, revealing a significant cut¬ 
off dependence in the 3D-SPGPE results, as shown in Eig. 6. 
This sensitivity is essentially the same mechanism behind the 
jump in Eig. 3 for A = 4: increasing the cutoff energy above 
will pick up the first transverse excited state in the c- 
field description, and so some care is required in selecting a 
3D energy cutoff that generates a consistent projection into the 
transverse ground state. Eor A = 8.5 a choice of cutoff mode 
population 4^^ = 0.8 or 4^^ = 1 selects an energy higher than 
required for the projection onto the transverse ground state, 
while 4ut - ^ gives a c-field density in close agreement with 
the ID-SPGPE. Eor A = 14 the behaviour is very similar, but 
now convergence to the ID-SPGPE regime is already appar¬ 
ent at 4ut - i ’ convergence in the tails of the distribution 
is seen for 4,,f = 2. 
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